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Two and a half centuries P. C. Archimedes li:i6/?-k;12) wrote 



a work entitled hE llS QUAE VEKUIJTUR IN lIUIvlIDC. He maintain- 



ed that every particlJ of a fluid mass,whbn in equilibrium, 



is equally pressed in uvery direction. The lav/s and j roper- 



ties of liquids were investigated by this ancient mathemati- 



cian sufficiently to unable him to aevise a h:drostatic 



mcthoa for aetermining the purity of the precious m-tals. He 



maae a further practical application of the results of his 



studies in this airection by the invention of the screw en- 



gine, or the Arcliiiiituean screw. 



The Alexanarian School is accreuiteu with the construc- 



tion cf a fe»v hydraulic machines such as the siphon ana the 



force pump. However fluid motion was probablj first stuaied 



by a Roman, Sextius Julius Front inus, inspector cf the j^ublic 



fountains at Rome in the reigns of Nerva ana Trajan about the 



close of the first century A. D. In his wori; cntillea ^^Di. 
AQUAEDUCTlEUn UREIS RCi..A£ CCMUiiNTARlUS, is I'ouna v uu scrip Lion 
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or the great aqueaucts of Ronie. Front inus construct ea fivc 



new aqueaucts,. ::uiKins in ail fourteen, l.e also at^scrjbes the 



methoas used in aetermining the amount of vater aischc-rrjea 



from ajutages and the methoas of aistributing the w^ater of 



the aqueaucts ana fountains. He observca that tlie amount of 



wrater aischargeu depoaas upon the height of the water in the 



reservoir abovt; the orifice as well as upon tl:e area of the 



orifice. 



Fifteen centuries later Castelli lloL.3 ) aavincea the 



theory tiiat the velovity of uischarge is } roj-ort ional to tiie 



heip:ht of water in the reservoir above the orifice. Torri- 



celli,a contemporary of CCcStelli,observini; that a small jet 



of water rusliing froai an ajutage nse'j lo l. ii-i^x.t rlmost e- 



qual to that of the water in thu reservoir , conduced that the 



velocity of discharge io equal to th-^ velociti which the 



water wcula acquire by falling freely through a aistimce the 



same as the heiglit of water above the orifice. As a result 



of this observation he obtained the correct relation between 



the velocity of discharge ana the heaa of water llG^o), 



Marlotte 1162C?-1G84) ,the author of a posthumous v.ork, 



entitiea TRAITE DU tICUV£M£NT D£S £AUX,maae extensive use of 



the theorem of Torricelli. He seems to have be^n the first 



to attempt to reconcile tlieory una exjeriment by attributing 



the retarac.tion to friction. The filiments sliaing along the 



surface of the pipe were supposea to be retaraea ana other 



filiments having a greater veioci^j' than those naer the sur- 



face were retaraed by rubbing against the slower ones. Tiie 



retaraation was supposea to be proporticncl to the ui stance 



from the axis of the pipe 



Guglielmini, a contemporary of Mariotte ,aevotea himself 



to the study of the motion of the water in rivers and csnals. 



Ke assumea that every particle in a verticle section moves te 
with a velocity equal to the velocity of aisclierge from an 



orifice at an equal depth below the surface of tiie water, and 



the coluji.n cf the fluia x.aving the area of the crifjcu as a 



bauiQ. Therefore this weight which is i-roport ional to th^ 



ea of the orilice multiplied by the height of fluia in the 



vessel must be proportional to the quantity of niotion engen- 

aered in the particles which escaje through the orifice each 

instant. Continuing tiiis analysis he finally t^rrives at thu 

theorem alreaay ostablisheu by Torri celli ^tt an entirely aif* 

furent n.ethoa. 



The investigation of the flov; of water in rivers attract- 



ed consiaerable attention in Itcly ,i.robably on account cf the 



extensive lanascapv:; garaening. pesiaes (Juglielmini, v/ho was 



insiectu^of rivers in the j-p lanese, Marquis Poleni aeserves 
mention in this connection. In IGJb he wrcio Di: MCTU a4q,UA£ 



MTXTC;and in 1718 another work concerning the flow of ifcater 



through orifices ana short tubes. 



Newton I164ii-1/L7) invest 1 gat ea the effects of friction 



ana viscosity in diminishing the velocit:v of running water. 







Ik 



6 



In book 2,§j,Gf the PRINCIPIA he offers the folioving hipoin- 



esis: The resistance which crises from the viscosij of a flu- 



ia, other things being equal, is i^ro^ort ioaal to the velocitj' 



with which the j articles separate from eacli other. Tiiis 

te 

may be saia to be the funuamental principal unaerli ing ell 



that part of hydroaynainics which deals with viscous fluias. 



The Vi.NA CCMTRaCTA as well as the oscillation of waves serins 



to have been consiaerea first by Newton, 



Daniel Fernoulli publisheu his hTDRCDYNMITC/^STVl DlC.VlR 
IPUS ET i\:CTirUS PLUIDCRU?4 CCMM£NTARn,in 1736, in which he 



bases his theory upon the suppositions that tiie surface of a 



fluid containea in a vessel which is buing emptieu by an ori- 



fice remains always horizontal ana that horizontal strata al- 



ways remain contiguous to each other, ana that the particles 



descena vertically with avelocity inversly projortional to 



the horizontal section of the reservoir. His principle was 



not acceptable to his contemporaries, consequently John Fei?- 
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noulli and Maclaurin oach atiemptea to so^ve the rroblem by 



inaependent methoas but without markeu success. 



Jean-le-Rond D'Alembert 11717-1783) aiaud by the aiscov- 



eries of Kuler 11707-1785) tooK the first great step in aeter- 



mining the general equations of motion of a purfiect fl':iia. 



When generalizing James Pernoulli's theory of pendulums he 



discover ua a principle of dynamics so simple and general that 



it reduced tlie laws of the motions of bodies to that of their 



equilibrium. Ke appliea this principle to the motion of flu- 



ids, and gave a specimen of its application at the end of his 



Dynamics in 1743. It was more fully developeu in his TRAITE 



DhS PLUlDilSjV/hich was publishea in 1744, where he resolves, in 
the most simple ana ellu^ant manner, ail the jroblems which 



relate to the equilibrium and motion of fluias. Ke makes use 



of the very same suppositions as Daniel Fernoulli, though his 



calculus is established in a very different manner. Ke con- 



siders, at rvery instant, the actual motion of a stratum &s com- 
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posea of a motion which it had in the preceeding instant and 



of a motion which it has lost. The laws of equilibriuii be- 



tween the motions lost furnish him with the equations which 



ug/i 



represent the motion of the fluid. Althcug}(the science of 



hydroaynamics had then made consiaerable progress, jet it was 



chiefly founaea on hypothesis. H remaineu a aesideratum to 



express by equations the motion of a particle of a fluia in 
any assignea direction. Tnese equations were founa by D'Aim- 



bert from two principles: 1st that a rectangular can8l,t£ken 



in a mass of fluia in equilibrium, is itself in equilibrium; 2a 



and that a portion of fluia, in passing from one place to a- 



nother, preserves the same volume when the fluiu is incompress- 



ible, or dilates itself accoraing to a given law when the flu- 



id is elastic. Kis verj ingeffiious methoa was jublishea in 



17b2,in his ASSAI SUR LA RESISTANCE DES PLUIDES. It wis 



brouglit to perfection in his OPUSCULAR MATKEKiATlQUES, and was 



adoptea by Euler." Philosophers fiad attempted in vain to ue- 



t ermine the laws of fluia motion independent of all hypothe- 



ses. However the method of flioxlons iroveu inaaequate to the 



purpose anu it was only after £uler haa conti'ibuied to sci- 



ence his calculus oi' partial differences that the object was 
reachea. D'Alembert first appliea the new calculus to the 



motion of Aater.and he anu Euler both succeeueu in obtaining 



equations of motion for a perfect fluid restrict ea by no par- 



ticular hypothesis. 



Chevalier Dubuat,a French engineer ,publishea in 1786 a very 



satisfactory theory of the motion of a fluia, basud upon the 

experiments 

experience of himself and othBrs. "Ke considered that if 



water werea perfect fluia and tlie ehannels in which it flows 



infinitely smooth, its motion would be continnually acceler- 



atea,like that of bodies descending on inclined planes. Put 



since the motion of rivers is not accelerated ,and soon ar- 



rives at a state of uniformity , it is eviaent that the viscos- 



ity of the watar ana the friction of the channel in which it 
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descends must equal the accelerating force. Dubuat there- 
fore assumes it as a proposition of funaamental importance 3 
when vater flows in any channel or bed, that the accelerating 
force which obliges it to move is equal to all the resistances 



v/hich it meets with, whether they arise from ti.e viscositj or 

from the friction of the bed." 

C 
In 177^ Abbe Bossut conaucted very uxt endive invest iga- 

t ions to aet ermine the amount of retaraation in pipes uue to 

adhesion and friction. Ke reachea the conclusion thc'.t with 



pipes of the same length the discharge is proportional to the 



dia*aeters,and when the aiameters are equal the discharge is 



inversly proportional to the square root of the lengths. 



The general equations of motion are so complex that Jit- 



tle could be done with them. This result ea in the numerous 



efforts which were maae by onii class of investigators such as 



Dubuat jPossut , and Prony to develop by experiment sim.pler form- 



ulae for the velocity of running water. Ilov/uver the more 
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theoretical mathematicians, among whom were Lagr ange, Laplace , 
and Poisson,enaeavo»ed to obtain solutions of the equations 
of D*Alembert anu thus solve th«j various problems of hydroay- 
namics. These men together with other eminent mathematicians 

both of Prance and Great Britain aevoted much time to the in- 

\ 
vest igat ion of waves ana tides. 

Navicr I178b-18b6) was probably the first to offer tqaa- 
tions cf motion for a viscous fluid (ISk^i:.), He bases his in- 
vestigation upon the following suppositions: -the fluid is com- 
posed of ultimate molecules, act ing upon each other by forces 
which, wiien the fluia is ct rest, are functions simply of the 
distance, but when the molecules receue from cr approach each 
other, arc i..oaified by this circunist : nce,so that two molecules 
repel each other less strongly when they are receding end 
more strongly when they are approaching each other than when 
at rest; and the alteration in the attraction or re^pulsicn is 
supposed to be proportional to the velocity with which the 
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inoleculeSySeparatea by a given aistance,are aiproaching ox* 



receding from each other: a second supposition involves the 



symetrical arrangement of the particles. Navier aealt only 



with an incompressible fluid ana arrivea at the following e- 
quations: X-1/^ lap/dx-A(aVcix^au/cly'+ciiVaz*) )=2>ia/^t 
and the other two can be #riten by making similar changes in 
the corresponding equations for a t^^-f^ct fluid. A is a con- 



stant aepending on the nature of the fluia. The other terms 



are the same as those usea in the following jageii. 



Poisson (1781-1840) aerived equations not only for an in- 



compressible fluia but also for an elastic fluia in which the 



change in density is small. He treatea the subject fro;.: the 



standpoint of an elastic solia, supposing the fluid to be con- 



tinually beginning to be displaced like an elastic solia, and 



continually re-arranging itself so as to make the presure e- 



qual in all airecticns as is the casu i^rith a fluia at resi. 



Kis equations are writtenthus 
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X- VP I ^V/ ctx-A ( au/ dxVau/ay+du/az ) -B a / dx ( au/ Ci.x+av/ ay+uv/dz ) ) 
=^u/c)t with the other two to corresjond. 



For an incompressible fluid these agree with those obteinea 



by Navier since In that case the expression for cubical ex- 



pansion disappears. It is to be noticea that both of the fore- 



going methods involve a consideration of the ultimate mole- 



cules of the fluid. 

Larre ae Saint -Venant 117T^-1886) first succeeaed in ob 



taining the equations independent of any consideration of the 
ultimate molecules. Ke attempts to connect the oblique pres- 
sures in aifferent airections about a point with the aifferen- 
tial coefficients dm/dx,diVdy 5:c,wiiich express the relative 
motion of the fluia particles in the immediate neighbcrhoca oS 
the point, by assuming the tangential forct; on a plane passing 
through the point to be in the direction of the principal slia- 
ing along that plane. Ke then employes theorems by Cauchy^ 
connecting the pressures in uoy aifferent directions in any 
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fluid, to show that the tangential force in any direction a- 



lons a plane is proportional to the sliaing along the saaie 



plane resolvea in thj same direction* Hu published this the- 



ory in 1843. 



In 1845 Q. G. Stokes derived the same equations as Poi- 



sson but found that A=3P. Ke basea his investigation upon th 



three principal hypotheses: IstThat tlie difference between the 



pressures on a plane in a given direction passing through c^ny 



point P of a fluid in motion and the pressure which woula ex- 



ist in all airections about P if the fluia in its neighborhooa 



were in a state of relative equilibrium aepenas only on the 



relative motion of the fluid immeaiately about P;that the rel- 



ative motion due to any motion of rotation may bo eliminatea 



without affecting the differe* ces of the pressures sbove mon- 



t ionea*. (Stokes) Ke then arrives at the conclusion thAt the 



stresses due to viscosity are functions of the rat( ii ol iiirtiin. 



After speculating as to the molecular constitution of the flue 
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iu he arrives at^ lh<^ hypothesis that these are linear funct- 
ions* This Jiis his second hjtothesis. Fcr gases he introau- 
ces a thira: "VVi.en g gas is expanaing equally in all uirections 
the stresses P,Q,,a:icl R are the same as if the fluia ^era fric- 

tionless " In his report to the Triuish Association in 

1040 CN RfiCuNT R..SEAKCHE3 IN HYDRCDYNAjMICS, Stokes clsims tna t 

the principal feature of his investigation consists in eximi- 

of a fluid 
nating froi.: the relative motion about^anj particular point 

the relative motion \Arhich corresponds to a certain motion of 
rot at ion, ana examining the nature of the relative motion which 
remains; ana that the metiioa employea does not necessarily re- 
quire tiie consideration of the ultiinato molecules. 

The secona assumption of qtokes is not altogether satis- 
factory since it rests upon the supposition that the volccitj 
is small. 



In lb61 C. il. l-iejer a^^rivea the orainary equations ii' 



F=0 whicii,accoraing to other invest i^jators prion to thio time, 
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at-e true only for an inconipressiblo fluia. He begjns liis in- 
vestigation aft^r some preliminary consiaerat ions with the 



rollottTing statement : *The internal friction of a fluia takes 



place between the different strata of the same fluid, end is 



proportional to the differential-coefficient of the velocity 



along the normal to the plane separating the strata." Pat 



this hypothesis explains only six of his nine initial uxprus- 



sions. For instance if we take the three initial expressions 



for the retarding forces parallel to Lhe x-axis 



**nluivclx)^y CLZ, 



-n(du/ay)dz dx. 



-nlaiVcLz)ax ay, 



The second and thira reaailj follow from the hi.fothesis but 



the first has no aefinite meaning, unless we extuna liis hjp^-^^^ " 



esis as W. M. Hicks nas interpr'-tea it in the rritish Assoe 



ciation report fcr 18bi so that it reaas-"The friction on a 
small plane in a given ai recti on in the plane is proportional 
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to the rate of variation perpendicular to the plane of the 



component of the velocity in the given airection, wliilst there 



is a normal part , proportional to tho rate of variation per- 



pendicular to the plane, of th« component i eri enalcular to 



this plane. 



The subjcict has been investigateu by Stefe?n (18^)*:;) anu 



later by Maxwell, M, Levy,M. Kleit2,and Putcher without aauing 



materially to what had been done already. 



However the work of Maxwell resultea in the aetermina- 



tionof the constant called the coefficient of viscositj . Fol- 



lowing Maxwell numerous experiments have been made for the 



purpose of determining the value of the coefficient of vis- 



cosity. Helmnoltz ,Piotro.»'Ski, Maxwell , Meyer , L'na Poissuuille 



may bo naiiied as having maue the most elaborate series cf ex- 



periments. 

Lamb, in his TREATISE CN FLUIDS^, aerives the equations of 
motion for a viscous fluia by a methoa basea on those of St. 
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Venant and Stokes. Basset , in 1888,publishea A TPuEATISiC CN 



KYDRCDYNAJvIlCS in which he uses Stokers 'methoa. Both authors 



accept the second hypothesis of Stokes so that we have no as- 



surance that their equations hola except in the case of slow 



motion. Both fully relized this and restricteu the applica- 



tion of their results accordingly* Cf all moLhoas so far 



that of Stokes haa recwivea the most general approval. 



In the following paper we have derivea the equations 



for the motion of a perfect fluid in order that we might ox- 



tend the same methods as far as possible in obtaining the e- 



quations for a visccucs liquid (the aiscussion is limitea to 
an incompcessible fluid, or liquid). The soconu part is ae- 



voted to the consideration of a viscous liquiu* The aiscus- 



sion is based upon the definition of the coefficient of vis- 



cosity as agreed upon by exp v^rimental phy;Bieisis, Then by 



comparing the motion of a viscous with that of a perfect li- 
quid we are enablea to give a aefinite uieaninff to all nine of 
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the initial expressions obtained by Meyer. To obtain the e- 



quations referred to eyllndrical ana siherlcal coordinates 



ve have not usea the method of transformation but have aeriv- 



ed them by analysis simfalar to that usea in uetermining the 



equations v^hen referred to rectangular coordine^tes. Fy tiiis 



jp.ethod a definite meaning attaches to each term in 'the re- 



sulting equations 
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PART I 



i?:QUATlCNS OP MOTICN CP A PLRPECT FLUID 



D 'Alouibert •& Principle: Since L*Aleri:burt 's princitls 



is useu so frequently in '^nex follows, it ma^ be wcjll lo say 



by way of explanation of this principle that its ust means 



the equatin;^ of 



Trnpresseci Force - Pressure - £xpruss8d Foi'ce to the forces 



lost througii friction. 



An la..al Fluia: V/iLh this as witii ctlie. irohloms it is 



easier to consider first the iaeal conaition. An iaeal fl lia 



is one in which there is no friction; the particles move bj 



ono another without causing any mutual ret araat ion; hence A'lih- 



out any loss of energy. 



Tlie Derivation ol the iiqu&tions of Wot i en of h Porluct 



Fluid: To aei'ive the equations of motion of c fci-fect fluia 



is a comparatively easy matter. Let A D (Pig.l) roii'osent 



an w'lement of such a fluia; x,y,ana z be the coorainntes of 



cl 




Fig. 1. 



the point A;A g,A F,ana A K l.e representee by ox, a; ,£na uz re 



sject ively;X?Y,anaZ ti:u^ imprussea force reaolvea larpllei to 



the x,;y,and z axe8;u,v,ana w me con^ponenLS of thu vtilocit^- 



in the same direct ions; anu finally let p represent ihe pres- 



sure at the point lx,y,z) or A ana rn liie .iiass of the elun:unt 



From Lhe conaitionts of Che ]:roblein 



u = f ix,: ,z,t ) 



ana it follows immeuiatwly that 



'^u/'2)t = a^v ^t 4- ula'j/ ux)+ vldiv'oi)+ vvlav^2)= fj, 



This last is an c^xpression for ihe effective force or ^xpr.sa- 



ea force in the direction of tho x-axis. In s aiMili^r no- 
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tation ?) v/?>t , cf f^ ana '^v/?>t , or 1"^ rej-r^sent the effective 
forces in thu uirection or the y tinui z-axes rusj uci i vel^ . 
The total impressed forces are iriX,mY,ana mZ;iht total ci'ectiv 
forces mf^ ,mf^ ,and mil ;the pressures on the foces A C,A £,c:na 
A F are p dy dz,p dy dx,and p dx dz respect ivelj ; on th^. fa- 
c-s opposite the pressures are p dy dz + (dp/ux)ax ay dz, 

p dy dx + (up/az )ax ay dz,ana p ax az + (ap/dy)ax ay az. 

of 
Therefore the excess pressures on the opposite faces aru 

t dp/ dx ) ax dy az , ( ap/ az ) ax dy dz , and ( ap/ aj. ) ax ay az . 

Since the hypothesis excludes thti possibility of any 

loss through friction,ve havi^ novr consiaerea vll thL? forces 

involved. equating separately the fcrcus acting in the ai- 

rection of each axis by D •/-lernLert ' s principle , v/t:; obtain the 

following equations: 

mX - mf> - lap/ ax) ax ay az = C 

mY - n^-v - lap/ay)dx dy dz ^ C 

mZ - mt\ - ldp/dz)dx dy dz = C 
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Substituting for m its Voluc|^ax a: az in v/hich /^reirc?- 
sents the aensity of the iluJa,ana then aiviuing b; ax a\ az, 
tJne resulting 'ijquaticns are: 

mP- iy^e - ap/dx = C 

Y/^ - n/^ - ap/a: = C 



Z/^ - i^C - ap/ CcZ = C 



or as usually written: 



X - 1/ /^ ap/ax =^u/'3t 



Y - \/e cp/uj ='^ v/:^t 



Z - i/Y^ up/az = l)ii/2\ 



So far we have obtainuu tiiree relations between tiie five 
quantities a,v,w,p,anar . Wh^^n the jquavions ere jscu in t.^ 
plication to an incompressible fluid. /^ is const ant ; but . wiien 



other fluias arti under conslaeration, a fourth equation is isui^ 
pllea by Boyle •s. Law, 

P = k/^ 



Continuity: Th- ccnaition cf continuili lurnishes a 



^'± 



fifth equation. If tlu fluid l)t cent inuous, then tho incruast.- 

Vw 

v/ithin a given Sj.ace aurinr; any given time /ust equal iliu ox- 
cess of inflow over the outflow across tiit bounuai-ies in the 



givjn time. 



The mass of fiuid containea by the element at any momt;nt 



IS 



pdx ay dz 



and its increase in tlit; time ui is 



(a P/dt)at ay ax az 



The flow ovef the surface ay az in time at is 



u P at dy az 



and the flow over the opposite face is 



u p at ay dz + u{\x (^ )/ ax at ax a; az, 



giving an increase i.\ the element of 



laluP)/ax)ai ax ay dz. 



Similar ilj tho increments in the airections of tiie i ana z- 



ax.^s are fo^na to be 



2b 



(a(v/^)/ay )at dx uy az 



anu 



(alw/^)/a2)dt ax ay az. 



Equating the \vhoiu increment to thv. pariiLl increments &na 



dividing by at ax ay oz v/e obtain 



d /^/at = atuP)/dx + atv(^)/ay + a(wrP)/az. 



Fut for an incompressible fluiq^, /^ being a coniittmt ,this be- 



comes 



du/ax + av/uy + dv/dz = 0. 



Cylindrical Coordinates: The equations of motion of a 



perfect fluia wJien referrea to cylinarical cooruinaLes raaj be 



aerivea from, those already obtained for rectangular coorai- 



nates by transformation. The following however is a more sat- 



isfactory mttiiod Jn so . uch as i t. gives a aefinite inter] re- 



tat ion to each terai appearing in the I'esulting equations. 



Moreover this methoa to be appliea an the aerivation of the 



equations v/hen referrea to cylinarical ana polar coorainaies 





k 
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forms the basis of several important princiiles involvea in 



the derivation of ihe equations of motion of a viscous liquiu 



referrea to th.; same systems of coordinates. 



The velocities u,v,and \ir,and the impressed forces X,Y, 



ana Z have thw airections of r, 6 ,ana 2 respectively. v is 



angular velocity. 



A 




fyFig.2 



X p r (3l6 dz dr = ilie impresses force in tlie airection of r 
(^u/St)Pr aB dz dr = efective forct; in the same direction 



due to an acceleration in the velocity u. The velocity v th- 



rough centrifugal force contributes to the effective force 



in this airection 



-V /^r a 6 dz dr. 



The velocly w can have no influence uj on the foi-ces in the 



L7 



direction of r. 



- 1 dp/dr ) r d ^ dz dr 



is the difference in pressure on the two faces C K ana D M. 
Equating by D *Alembert 's principle we have 
XP r d& dz dr - (dp/ar)r a 6 dz dr - tT)u/^t)/^ r a© uz dbi 
+ v^^r^dfi dz dr = 
Dividing by p r a dz dr a.n{x transposing, this equation be- 



comes 



X -Jj^up/dr = '^u/^t - vr. 
Similar ily for the forces acting in the airection of ^ 

Y p r a 6 az ar + impressea fcrco, 
rOv/^t)/^ r a© dz dr = effective force con- 
tributed by acceleration Jn v. 

ll/r)(dp/de )r dfi dz dr = difference in the 
pressures on the opposite faces. But there is still another 
resulting forire acting in this direction aue to the fact that 
while a particle is being carriea in the direction 6 it is 
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also being carried in the cirt-cticn of r bj a veiocitj 



u; consequently at the end of the time dL it will have fallen 



behind by an amount equal to 



(r + dr)a^ -r d6 or ar dQj 



which may be written 

r d^ dr 



r. 



Prom, physics 



s  1/2 ft* 



f= (2s)/ f^ 



where f = acceleration, ana s the space aescribed and t the 



time. Subst itut ing 



f = (2 ar ae )/at^ = 2uv 



which Is the force due to the fall. 



Multiplying by (^r a© dr dz we get for the total of this 



force 



^uv(^rae drdz. 



The velocity along z as before contributors no component to 



2d 



the forces in this direction. 



Equating the forces founa ana simplifjing w'e get 



Y - n/P)[l/r){up/de) = rl^v/St) + .:uv 



Neither u nor v can contribute a component to the fore® 



acting in tlie direction of the z-axis, therefore the equation 



of these forces will be id.^ntical with the corresponding e- 



quation for rectangular coordinates. 



Collecting the results just obtainea ¥/e have 



X - ll//^)(ap/ar) = Zlu/at - vr, 



Y - ll/r)ll/p)laF/ae) = rOv/pt)-^ iiuv, 



Z - H//^) lup/az) -^w/dt. 



The equation of continuity is 



d(ur)/dr -•- rldv/aa) + r law/ uz) = 



It is unnecessary to aerive this last equation since the meth- 



od by which it is obtainea is iaentical with that used to ob- 



tain the similar equation referrea to rectangular coordinates 
Spherical Coorainates: Angle WCP iFig.b.) ^6 ana 
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ig.o. 



PCD = Q> . Let u,v,and nr represent the velocities in the ui- 
rections r,0,ana (D rosj ect i vuly ;X,Y, anci Z ^v^ symbols for 



the corresponaing impressed forces. v ana v/ are angular velo- 



cities. The mass of the element P .M is 



/^ar rd 6 rsinefi a(D = m. 



Xm = impressed force 



in the airection of r; 



(3u/^t^) J^ = effective force 



due to the acceleration in u;ihe cjntrif^igal force 



-r 



v'm 
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results from the velocity v;tlie centrifugal force contribuiea 



by w is 



-r sinu B w^m 



in the direction of D and this gives the component 

-r sinci^fi v/ m 



along the radius vector r; 



(ap/ar)ar ra^ r sine 6 aO = 



the aifference in pressure on the two opposite faces. 



Equating and transforming as before 



X - ll/P)lap/dr = '^u/2>t - r v^ -r sinc^ w^ 



With ruferenctf to the airection of ve have 



Y m = impressea force, 



rOv/^t) m = aifferencu in pressure 



on the opposite faces. 

Since as a particle moves from e to p it is carried by u tlie ' 



distance e R it falls behind by an amount equal to ar a6 



vrhich when expressed in terms of force is equivalent to 
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2 u V 



or multiplying by mass 



2 u V m 



represents the component or effective force coniributea by 



the combined action of u and v. w cannot exert an effect ivo* 



force such as this since the particle always moves in the 



same parallel, but it does exert a centrifugal force in the 
direction of ^ ,or r\ather there is an effective component of 
the centrifugal force acting in the direction of 6 due to w. 
y/e already know (see page 31) that the centrifugal fcrce aue 



to w alon^ D is 



-r sinu 6 w. 



Since the tangent to e F at e makes in angle G witi. C D,ihe 



component of 



- r sine © v/* 



along the tangent is 



-r sine 6 cos 6- w. 



2>b 



Multiplying by mass we have 



-r sine e cos £? w m. 



Equating and simplifying, the I'ssulting equation is 

Y - ll/(^) (l/r)(ai/d^) = rd^v/dt) -hk^uv-r sin.P cos a wf 
The forces along the O coordinate are 

Z m = ih'pressed force 
slne^ ('3w/^t)r m = effuctivv^ force 
due to the acceleration of tiie angular velocity w; 

ap/tpr sinefi d^ ) m = difference in the 
pressure on the opposite facos. 

Since by the combinea action of u ana w a j-Q-^ticit ic^lls 



behind 



dr sine 6 aO) 



vhichjConverteu into the terms of force, ^ives 



2 sine 6 a u. 



Therefore we obtain a total force cpntribui.a b> the combin- 



ea action of these two velocities equal to 



o4 



2 sine © w u m; 



V also carries the particle through the aistencc e F anu sinoi 



the line P K is shorter than e K there will bu a resulting 



effective force aue to the combineu action of v and w 



e H = r sined a ^ 



P K = r sinc(0 + a^) a <I> |see note \3(,) 



= r dfl) (sinefi cos ^^ + cos^ sin^^^). 



Substituting 1 for cos a^ anu a© for sine a^ 



P K = r a© since + r u« cosfi d* 



P IC - e I- = r a«P cos ^ a^ 



or the aistance tiirou.^h v/liich tiie particle ma; be saia to faJL 



unaer thtr combineu action of ^ anu v. This is uquivclent to 



a force 



2 r cos ^ w V , 



or for the mass of the elem^. nt 



ki r cosd w V m. 
L 

Py equating anu simplifying we obtain the tl;iru t;4uation 
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Z - dp/(pr sinu# d9) = r sintBOw/jt) k:, sino^ u .v -f-Lrvw cos^f 



Equation of continuity: Thu increase? in thci t.lc\..onu 



fci* the time dt is 



,2 ^. 



(ap/at )r^ sineO ci(J ar ufi; 



thfc oi fferencu betwr^en the inflow and the outflow in the ai- 



rbction of r is 



(alarV )/ar)sinue a (Pdr ud; 



in the airecticn of 6 



(d(rv /^ sined )/ Ir a6))r^d0d&(xr; 



anu in the airecticn of (D 



(air sineew^)/lr sincf©a*))r sin^fi uiPa^ar. 



Equating ana aividing by r afl? a^ar 



laf^/at)sint © =(alur^p )/lr ar))sinb^+ dlrr /^ sino^)/lr dO) + 



la(r sine©w/^)/lr sine ^ aO) )sine ^ 



or 



(d/yat)sine ^ =(a(ur VVl^^^dr ) )sino 6 -«- dlv /^ sine a )/a e + 



(dlw p)/d 9 )sin«e 






or \t P is a constant 

(d(ur^)/(r^ar)j8ine^ + a(v sinee)/G.e +(aw/a © )bine^ = 0* 

Collecting the equations obtained ^e have 
X - dp/(/^dr) =i>u/^t - r v^- r sine'^w^, 
Y - ap/(/^r ae) « r(^v/^t) + 2uv - r sinse cosa «r^, 
Z - dp/(/^r sine 6 dtf)) = r sine© {5«r/3t) ■«- ;ii sine^uw -•- k;r cos^ v«r, 
(d(ur^)/ (r^dr ) )sine + ia(v sinee)/de + taw/a^)sine6= €• 



NOTE; Cn page 34 v/e say P K = r sine(fi+ u&) a(P, It 



is to be noticed that d© is taken in the nec^ativt; sunse. 
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PART U, 



THE EQUATIONS OP MOTION OP A VISCCUS 



LIQUID 



Having determined the equations of motion of a j-erfect 



fluid which is incompressible cr compressible when referreu 



to the principdLi systems of axis, we will now proceea to the 



consideration of an incompressible £±ut viscous fluid which 
we shall call a liquid* In part I. we have saia little about 



a compressible fluia in the aerivation of the equations but 



a littli consiaeration will show that the same form foB the 



three principal equations holas true for both kinds of fluias. 



However in what follows we shall have occasion to refer to 



the foregoing equations as applying to a perfect liquiu|gnly« 
The problem we are about to consider involve^ the forces al- 



ready dealt with in the case of a perfect liquia together 



with the forces lost throuf^h internal friction. 
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The problem is to aet ermine the relations existing be- 



tween the forces lost and the velocities. 



Coefficient of Vfftscosity: It is necessary at this point 



toa^fine a quantity callea the 



cotjfficient of viscosity. 




Pig. 47 



If A £ D R represents the initial position of an element sub- 



jected to a shearing force exertea parallel to the planes of 



the strata e.g. plane A T R B such that the element takes the 



form. A M G R, the coefficient of viscosity is aef ined as the 



ratio of the shearing force to the shear per unit of time. 






Putting this aefinition in the form of an equation we have 

(shearing f or cc )/( shear per unit of timu) = /^ 



in whichy// is the coefficient of viscosity: or shearing force 



-ylf times shear per unit of tim.e. 
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Let u be the velocity in the direction of A B;and the point A 
be located by x,y,z;ana the sides of the element be dx dy dz. 



Then 



u = f(x,y,z) 



du/dy = f*(x,y,z) = 



rate of normal displacement along y. This force acts upon 



the surface of the element dx dz. Hence th'e shearing forcb 



on this surface is 



P =y^ (du/ay)dx dz 



The force on the opposite face is 



F' = /t/ [du/dy) (XX dz + y^ (ali/dy^dx dz ay. 
The difference of the forces of the opposite f^ces of the el- 



ement considered is the retarding force due to the viscosity 



of the liquia. 



Derivation of the Equations: Sufficient has now been 



saiu to proceea to the derivation of the equations. ?/e shall 



take up this problem in much the same «ray thai wc aia the 



k 
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problem of a perfect liquia# But this ne« eibment of viscos- 



ity must be taken into account. It is tiviaent that our equa- 



tions, when derived, must be tho same as those alruauy obtained 



for a perfect liquia with terms attachea to account for the 
retardation aue td viscosity. For if in our neu cquations/V 



be equatea to zero the equations for a perfect liquia must re- 



main. 



Forces to be Equated: The forces to be equated are 



Impressed force. 



Pressure, 



Effective force, ana 



Retarding forcu,or force lost, 



which give when equatea 



Impressed f orce-Pressure^iUff ect ive= -Pjtarding force 



or in the more usual form 



Impressea f urc«-Prcs8urt=Erf ect ive force-Ret araing force. 



It is to be noticed that the sense in which tlie retaraing 
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force ia taken is opposite that in which the impressed force 



J8 taken. 



Rectangular Coorainates: 

1 



/ 



/o 




X 



/ 



/ 



l 



Pig*5. 



Take an element p b,the coordinates of the point A buing x,:^, 



and z ana the uim^nsions of the element being ax ui oz. 



Suppose the velocity uniform in the direction of x. 



Then any parallelopipea such as a b woula m.ove as a solia;bui 



the element is composed of a great many sjuch parallelcpipeuti. 



If the velocity in the airecticn of x varit;s as we pass from 



A to K,the layers in slipping by uech other will proauce 



mutual retaraation. The rate of change in u along the normal 
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to its direction is 



aiv^ay. 



From the definition of viscosity it follows that the shearing 
force upon the face s b is 

^x^ (du/ciy )ax dz 



and that upon the opposite face is 



r?, > ...ty 



^^(du/ay)dx dz + ^ { au/ ay) uxi. ay dz 



indicating a loss of force within the element equal to 



^ (au/ay)dx dz dy. 



In a similar manner we obtain the retardation upon the face 



ax ay to be 



,.^(du/dz)dx dy 



and on the opposite face to be 



^(au/az)ux ay +^(du/a:?)dx dy dz 



showing a loss within the element of 



^ ( d V d5?) dx dy dz 



All of the loss within the element is now accounted for pro- 
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vided the velocity in the direct ion or x j s nnilorLi. It re- 

vtlocity 
mains to examine the effect of an accelerating £^f^e^ alon^ x. 






To simplify the problem suppose the motion two uimeny^ional. 
In this case it is only necessary to consiaur a section such 



as A C in Pig. 6 




Pig. 6. 



•■ 6 i r 

Suppose unear certain conditions a certain motion exists at B 

will 

Tf the liquia is perfect a certain amo:;nt^t5* flow into ihe 

parallelogF^- A C across the line A E ujie to the acceleration 
of the velocity u. £ D represents the whole increase of in- 
flow at £ over that at A. 



£D = AB-DC. 



(A P)/ (D C) = uVu 



where u is the velocity across A B anu u* that across D C. 
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By composition 



lA E - D C)/A B = lu*-u)/u' 



which i-epresents the uistance of inflcv/ at x. when a unii i^ngih 



of A B is considered. The liuie of inflow is dt honce 



lA E - D C)/(A B at) =lu* -u)/ (u*al) 



gives the rate of aisivlacement in the aire ct ion of x. £xpanu- 



mg 

ll/at)((u'-u)/uM = ll/at)ll - u/^^M ^ 
(l/at)(l - iV(u-l-au)) =ll/at)ll-(l - u V^ ■»■ au^u'^ dui^u^+ uu^a^eic;) 
= ll/dt ) (du/^- auj/u^ +au/ u-^- uu/^ o^ ut c. ) 

Since tiiis sories canverrjes very rapiul^ udiless the incrwaae 
in motion is very great ail of the ttjri.is exceit the first ixiaj 
be droppea* This gives 

a\V l^ ^t ) = au/ux 
The retarding force on the face uj . az is 

yc/ I au/ ax ) a j az 
ana on thv3 opposite face 
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yH (au/cix) ay oz 4- y^ ( dS/ ux^ ay az ax, 



showing a losa for the eleraent of 



ytf (du/dx') dy dz dx 



when the liquia is viscous. 



The retardation in the airection of khe x-axis must 



equal the retardation of inflow, for there woula be aiscontl- 



nuity if it were otherwise. 



It is necessary to extena this principle to the case where 



the inflow is across two faces instead of one. 

' — \ 




ab . aa- hi . ij = 



total areal increase of inflow. 
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( ab au)/lhi ij ) = uVu . 



By compoaition 



(ab ad - hi lj)/lab aa at ) » (u* -u)/ (u'at ) 



equals rate of inflo^y. But 



lu*-u)/lcLt u') = ll/at)(l - u/(u + uu)). 



Expanaing this ana neglect ing, as before all terms except the 



first we get 



aiVcix 



for the rate of inflow. This shows that the lost forces ob- 



tained for two aimentional flow hola in the case when the 



flow is over two faces. This coula have beun concluaea with- 



out the above proof Tor only so much can flow in ana it is im- 



material whether it all flows in the air^ction of z cr partly 



in the direction of z and-p^rtHrly vn— the—direct^rjn of -z ana 



partly in the direction of y. 



The forces in the airection of the y ana z axis can be 



) ^ 



obtained in the same way thtst tiiose for the airection of x, 
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It is to be noticea that a force element; of the secona 



order is neglect ea if , for example, ih*-re is a norinal increase 



in the velocity Ujfor in that case the inflov/ will be accel- 



erated. 



Equatinf^ the forces found ana writing the two remaining 



symetrical equations we have 



X - dp/tp dx) = aVSt -^( au/dx^-»- ^V c.y' + au/az) 
Y - ap/l p dy) = ^v/^t - ;^ dv/dx^+ dV/a/ + ov/dz^) 



Z - dp/ 1 pdz) =7>w/3t - ^( d5f/^+ dw/dy'+ aw/ti^) 
Cylindrical Coordinates: 



h 




i^ 



Pig. 8. 



In treating the equations for cylindrical ana apiierical 



coorainates it is founa tliat the totel rate of aispiacement 



as for example on the face C K along the normal is not ex- 
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pressed by dia/cir for if the velocity u is constant there liiust 



be an inflow equivalent to u/r. Consequently in the follow- 



ing treatment account is taken of this anu then the element 



considerea as uniform in size i.e.r d a clz is a constant so 



long as we are consiaering tht. acceleration in u along r. 



With cj^lindrical coorainates for the tv/o analogous cases, 



the acceleration of v along 6 ana w along z,the opposite faces 



are equal in size; but in spherical coorainates in t^o of the 



three cases, u along r anu v along ^, the opposite faces are 



different. These turo cases will le aealt with as the con-e- 



spcnuing case in cylindrical coordinates 



One more consiau^rat ion is necessary. It has been sjtiown 



that V contributes an acceleration in the direction of r equal 



to 



-r V. 



If the change in v aue to a change ut in the time is av 



then the force 
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5 

r V 



gives arate of aispiacement at the surface or dz in the ai^ 
rection of 6 equal to 



-k. V 



!•©• angular displace nenl . In a si/:ilar manner we obtain 



two rates of displacement from 



ki uv 



which contribute to the forces acting in the airection of ft 



VIZ. 



2 u/r 



normal to the lace dr uz and 



2 v/r 



normal to the f ac^ r d e dz. 



Now we may write the rates of displact^ment . The total 
rates normal to the thwee races that contribute to the retaro^- 
ing forces in the airection r are 

(1) uu/ur + u/r normal to r aOaz, 
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(2) diV(r d^ ) - 2v normal to ar dz 



(3| uu/cLZ . normal to r ao ar. 



The corresponding rates for the direction & are 



(4) dv/dr + 2 v/r normal to r a^ dz. 



(b) av/de + ki u/r normal to ar dz. 



(o) dv/dz nonnal to r d^ar. 



And for the direction z we have 



(7) dv//dr normal to r ad uz, 



18) aw/(r da) normal to ar az, 



(y) dw/dz normal to r a ^ ar. 



Remembring that the surfaces over which rates (l),(b),aiil 
(y) act are to be consiaerea constant --(two of them are actu- 



ally constant )-- ana multiplying the rates by their respective 



surfaces ana f inaing the aifference in displacement for each 



opposite pair of faces wy have the following nine aifferences 



corresponaing to the nine initial oxprc^ssions on pages ^y ana 
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>1 



(1) 



(2) 



(3) 



14) 



(b) 



16) 



(7) 



18) 



19) 



du/clr^+ <Xu/ (r dr) - u/r*)r a az ur, 

au/l^'o-^^) - (2av)/ (r u#))r d 8 ar dz, 

dti/az')r d fl ar az, 

dV/ar^+ (3 av)/(r dr))r d a az or, 

dv/lJ^ ci^^) + (i:. uu)/(r^a^))r d ^ dz ur, 

uV/^)^ a 6 dz ar, 

dw/ar^-f uw/(r dr))r d ^ az ar, 

d^/(r^d© ))r a ^ az dr. 



2 

aw/az^)r d p az ar. 



The differences obtalhed from (^) anu (6) are rates in 



change in angalar aisplacement as comparea witii linear dis- 



tances ana consequently must be muliipliea by the raaius vec- 



tor before being equatea# 



J.ultiplying these expressions by the coefficient of vis- 



cosity, anu then aauing them to the corres|.onuing equctions 



for a perfect liquia wa have 



X - ap/l/^dr) = ^u/H - r v 
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/The letter/^ used in these equations equals ^ the coefficient 
of viscosity aivided. by the density. 



Spherical Coordinates: 




Flg.;i. 



(1) -r V contributes rate riiv rcrMal to r sineeuOar 



(2) -r sine*(9w^ 



• . " -ii sine6vv • r a^ar, 



(o) 2uv !! 2u/r./ r sinL 6 a^ ar, 



(4) 2uv 2v/r 



r a ^ r u(P sine^-, 



(5) -r sine e^ cos ^w^...», 2 cos e v/ r d 5 ur. 



r 




(6) 2 sine $ \'i u (2sine t)/(r sinea), . . r sines c. r d^, 

(7) 2 sine e V u (^.sine ^ u)/ (r sined^),.. r u^ur, 

(8) 2r cos d \ w (id CO 36? w-)/binu^ r sino 9 dO uv ^ 

[^) 2r cos a V w (L cos a v)/sine0 r ci6>dr. 



In taking account of the acceleration of u along r ana 



of V along 9 it must be reinei.ibured that if everj thing is con- 



stant there is an inflow in the first case equivalent to 



21^5^ 



and in the secona casj equivalent to 



cot V. 



In dealing with si:herical coorainatas we shall consiaer the 



element of uniform size in the same caseis as in cylinarical 



coorainates. The rates of aisplacement are 



ll) du/ar + 2u/ r on face r sine^d<Pd^, 



(2) uu/ (r ae) - 2v r bine© a<Pdr, 

-2sine w 

(3) aj/ Ir sine e od) r d^dr, 



(4 ) dv/ar + 2v/r r^ sinJ uQ a^, 



(5) clv/(i^,+ ^m/r -f V cot 6 on face r sine & uW ar, 



(G) av/ (sine £^ cl<p) - 2 cos ^ w r d^or, 



(7) dwr/dr + 2¥/r r^sinod^ u*^ a^, 



(8) dw/d5 -f ^c.ot^'W r sincd) d<?ar, 



i9) aw/a(f + 2u/3^ + 2cot^ V r d 5^ dr. 



Remembring that the surfaces jpon uhich (l),lij),and t*:^) 



£ct are to be considerea constant we obtain for the uifferen- 



ces in displacement 

ID taVcLr^+ 2du/li^ cir) - 2u/j^**)i'^ sine© a(f a^ur, 

(2) (aVlr^ciO + du/lr'dC'jcot^ •2dv/ (r ap ) -2v /r coi^r sine<^ 

(31 (du/lt*^ sino'& dtfM - 2dw/(r al? ) )r^ 3ine ^ a (^ ar a<?, 

(4) (av/dr*+ 4av/(r dr) +2v/r^)r^u^ aO ar sino-d, 

r^ sino' 6 uo ur u^, 

(5) (av/ (r a»*) -f 2au/ (r^ a^) + av/ (r u6;) cci^ - v/ (r sinu^)) 

(6) (av/(r sine^a a^P^) - 2cos ^ d\v/(r sinud a(P))r'sin-^ afl^ ur a^ , 
(8) (uw/ ur + 4 u»r/(r dr) + 2;Y/r^)r' sin^e aip.uv ae, 

(8) (aw/(^ clO -f 3cot^av//(r d^) -2w/r )r' sinj <^ d ^^^ dr d^, 

r^sine^d^ara^^ 
(y)(aw/(r sinee a<p') -f 2 cLV ^^^ ®^^^'^^''" 2cct^uv/(r sin^.^ u0)) 






The expressions deriveu from (4), (7), and (8) : lust be mul- 



tiplied by the proper radii before equating. 



Finally we have by annexing these quant it ic s, mult ipliea 



by the co^efficient of viscosity, to the cort-jsponuiag ^Cita- 



tions for a perfect liquia 

X - dp/(Pdr) ='3u/^t - v^r - r sine^fiw^- /^(du/dr'-»- 2au/(r ar) 
- ^.TVr^+ du/(v^a^^) + cot^ u\Ji/[r'aB) - k.uv/(r d6) - cote k.v/r 
4- dV(r'sine'e uifi) - 2aw/(r a^:^) ) , 



y - dp/(pr dd) = r Djr/dt -hiouv - r sined cose v^-/^(r av/dr -h 
4 dv/dr + 2v/r + dV/(^ cL^/') + l dv/ (r^ uB) -f (dv/ (r d<9))cot5 - 
v/(r sine^o) + 6y/{r sinea d?') - ^^ cosf? dw/(r sine^ a<P)), 
Z - dp/(psine0 u9) = r sine ^ 'dvr/^t +2 sine 9 uw -f 2r cog^ v w 

i. 1 9 2 

- /^ ( r sinea dw/dr -•-48ine© dw/dr 4-2 sine^ w/r +sine^ dw/(r u© ) 
+ 3cos e dw/(r de) - 2w sinea/r + aw/(^ sinu^ u(P) + 2aiV (r^ sire e df) 
+2cotfidv/(r sinee a(j)) 



Colorauo Si rings, May 26,lijjG--commonceu J^in. ,1(J;jI>, 
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NOTE A. 



For convenient reference the various equations are col- 



letted from the different parts of the foregoing paper; anu a 



tabulated analysis of the terms appearing in the equations is 



aaded. Special attention is given to the terms in the equa- 



tions referred to cylindrical and spherical coordinates. 



EQUATIONS FOR RECTANGULAR COORDINATES 



For a perfect liquid; 



X - dp/(Pdx) = )u/<)t, 



Y - dp/(Pay) = ^v/3t. 



Z - dp/(Pdz) = "Jwr/O t. 



For a viscous liquid; 



X - dp/ (Cox) = ^yj/dt - (du/dx + diV'ay*+ au/az ), 
Y - dp/(^dy) = "^v/^t - (dv/ax'+ ait/di'^-^ dV/dz* ) , 
Z - dp/(/^az) = ^w/^t - (dw/ax^-f aV/uy* + dw/dz'). 



Equation of continuity; 



du/dx + av/ay + av/az = . 



'^•i 



EQUATIONS FOR CYLINDRICAL COORDINATES 
For a perfect liquid; 
X - up/(f'ttr) = d\x/d\. - v^r, 

Y - dp/((^r (Id)rr V'<^t + 2uv, 
Z - dp/l^dz) = dn/h. 

For a viscous liquid; 

X - dp/l(^dr) = iiV^t-rv^ - ^(dVctr^ + au/' {r dr) - u/r*+ aij^/ ( i*' a 6*) 

r 

- 2 dv/lr (xd) + du/az)/ 

Y - dp/C/'r dS) = r^v/«)t +2uv - :^(rdV/ar* + 3av/ar + i^/ {v (i&) 

+ 2au/(r*a^) + r uv/dz' ) , 
Z - dp/(^az) = ^w/./t - ;:lldif/dr + aw/(r dr) + dvi/ [v^d'b) +• cilr/az'^ ). 
Equation of continuity; 
d(ur)/ar -f- r av/a6 + r aw/az = C/ 

EQUATIONS FOR SPIIEEICAL CCCEDINATES 
For a perfect liquid; 
X - ap/l/^dr) = £)i3/e)t - rv^- r i^ind wr^, 

Y - dp/l/'r a^) = v^/^t + 2uv - r sinfi- cos^sr* 
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Z - dp/(^r 8ind a<c) = r sin*^w/^t + 2 sind u w + 2r cos<? v w. 



For a viscous liquid; 

X - dp/l/'ar) = au/^t - v^r - r siAsv' r ^oh/^^^^ +'^ciu/ (r ur) - 
2u/r + du/ ( J^"^ cl6 J + cote aVl^^ ci6*) - 2av/lr a/9) - cot^ 2v/r+ 
du/l^^«ina dO?) - 2 aw/(r d^)), 

Y - dp/(/^r a©) = rhr/^t -f 2uv - r sin& cosew"^- ^(r dv/ct^^^+ 
4 dv/dr + 2v/r -h av/(r da^ + 2au/(r^d5^) + cot^ dv/lr a») - 
v/tr sinV) + dv/lr sin^ d<i^ - 2 cos9dw/lr sindd(^)], 
Z - ctp/l/^r ain9 ^(p) = r sin&2>w/^t + 2sin&uw -h .:r cocevw - 
^(r sing dw/dr + 4 ainadw/dr -f;. sin^i«r/r -h sln^dw/l^^ O-^) + 



i 



3ccseciw/lr d») - 2vir sin^/r 4- dV(r «in^ di/) + La\i/ [r^ sin/? acp) 

+ 2 cote- dv/lr ain&d^)). 

i':quat ion of cont inuit j ; 

a{ur^)/ Ir^dr ) sin6+ dlv sin&)/a©+ sinSavif/a^* 0. 



AN ANALYSIS CP TKi£ EqUATKNS 



Rectangular coordinates: 



The terms founa in tlie equations wjcien r(:;ferrea to rec- 



I 




t>9 



tangular coorainates need but littlt explanation. The terms 



dhi/dx\du/dY\a.nd uSi/az* 



are due to the variation in the velocity u in the directions 
x,y/and z respectively. And the terms found in the other b- 



quations are aue to a similar variation in v and w. 



Cylindrical coorainates: 



In equations for a viscous liquia 



In the airectio'^. of x-axis 



-r v'^ results from centrifugal force due to v, 



du/ar re suits from uu/ar , 

du/U' dr) - u/r results from displacement u/r^ 

Isee I age -SS) 
au/tr^a^ results from du/a^, 



-^ dv/tr d^) results froi-i centrifugal force -rv^ 



.2- 



du/az results from div'az; 



in the direction of y-axia 



2uv results from combinea actionccf v and u (se;.^ p Xf^ 



GO 



r dv/dr'+ dv/ar results from av/ar, 



I 



2 dv/ dr • 2 u v , 



dv/(^ cte^ av/de, 

2 du/(^^ d^) ^ u V, 



r dv/ dz ^ • • • • dv/dz ; 



in the direction of the z-axis. 



The equation Tor the forces actinj; in this airecticn is 



exactly the same as that for rectangular coorainates. 



Spherical cooBdinstes: 



In the airection of the x-axis 



-v^r reijulta from centrifu.^al force uuc to v. 



^ 2 

r s in ^ v/ • • • w , 



4 z 

au/ar tt^sulta from du/ar, 



2 du/l^ cir) - iiu/rf k.iVr (see i agu ^J) , 

du/(^^ciy) + cot e du. lr^d»). . . .au/a^, 

-:iav/(r d6') - cot# Lv/r centrifugal force-rv, 

au/(r^8in^ d^^ au/uf? , 
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k:a7//(r a(p) • • • .centi'ifucal force »i^ 



in the direct icn of the y-axis 



L. u V ccmbinec action of Ujcv, 

aue to nr, 
p srjn0 cose w^ • . • . .centrifugal force. 



^ 2 

r dv/cir + *:dv/cir • • • dv/ ar , 



2cLV/ar + ^v/r •••••••••• k, a v, 



2 » 

av/(r d6) av/d/9, 



k:, au/(r^a^) L u v. 



cot<9 av/lr d^) - v/lr sin<9)..,v cot^ (see pact^^^j-), 

uv/ (^ sin 6 a(/) .CaV/u^, 

-r sin^ cosdw"^ ; 
-k: cosdav7/(r sin<9 aC^). centrifu{^al force 



ill tiie direction of the z-axis 



26ind«r ii aue to cornbinea action of ujc*'lfie'>j i ./7 ) , 



k,r cos 6 V w v&v/lsco i^.J^), 



• 



« • .« 



r 8ineaw/ar*+ k.sirie dv//dr. . . .ciw/cr , 



) 



^sinft aw/dr + 2sin^ w/^^ >k:sin© w u, 



sand aw/lr a^) + co8 0aw/(r ad) . • • .aw/ud, 



2co3 t>^ aw/ (r aO) - sin rf 2w/r 2r cos^ v isf, 
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aw/(r sini^dt) ) uw/a(p. 



du/lr^ sine da?) ^ sino y/ u, 



2 cot^ av/ (r nine d<y) 2 r coa^^ v w. 



63 



NOTE B 



Advantages in the preceeuing methou of uerivation. 



The method of deriving the equations for a viscous li- 



quia 'AThich has been received most favorably is that of j^tokes. 



Foth Lamb ana passet have used his methoa,iJ:ie first in i art 



and the secona entirelj^^in their works on 1;^ aroaynainics, 

Eut StOKeti method is basod upon two ii^ poiheses. (so<- In- 
troduction). The second hypothesis has not been receivea 
very favorably because it involves the assuni^ion that the 
motion is slow. 



Tlie methoa which 7 liave used is basea upon the uefinition of 



the coefficient of viscosity; aru involves no lijpothesis 



save that of continuity of motion » 



The aerivaticn of the equations is much si.npler bj my 



mwthoa than by former m.ethoas unless an excejtion is maae of 



the methoa of 0. E. Mey^^r. The equations of Mejer are b- sea 



upon an assumption which aoes not satisfactorily exj lain part 



d4 



of Ihe terms (see £ntroduci ion) . 



The development of the equations of mot ion bj analjsis 



for cylinarical ana spherical coordinates enables us to give 



a definite ..eaning to each term in tne equations and may aia 



Materially in tneir solution. 



The iiiotiiod definitely aet ermines whon the equations hoiu 



ana when they do not as is shov/n in NCTJi C. 



The principal feature of this paper is tiie relation es- 



tablished betwoen the motion of a perfect ana a viscous li- 



quialseo pages 37 to 47). 
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NCT£ C 



ThS first apprdoximat ion which was maae in ihis paper is 



founa on pages 44 and 46 Jirhere all terms cxcejl the first are 



neglected in the expression 



(l/dt)(diVu - div^ + uu/u - du/u etc) 



This approximation is auiRij.>^.ible unless the accc;lerat ion in 



the velocity is very great which is not usually the case. 



In the seconu place v/e suppose the rate of inflow over 



the face in the direction porpenaicular to ihe velocity con- 



siderea to be uniform along tiie normc.;l to the face. It is 



easy to conceive that the inflo*/ might be acceleretea. Put - 



if it were th*j acceleration ne^lecteu woulu be of the seconu 



oraer and migiit reasonably be neglectea. 



If our conclusions are true the equations hola not onl^ 



when the motion is slow but in all cases L-xce] i where the 



the motion is i^apialy accelurateu in whicii catit. v/e coulo not 

so meny 
neglect t-tie terms in tiie equation above. 
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NCT£ D 



Cn the motion of a liquia in which a circular cyiinuer cl 



infinite length is slowly rotatinc* 



Whatever motion takes place in the liquia is aue to in- 



ternal friction rna is two aiinent ional, i.e. in the jlane 



perpendicular to the axis of the cylinder. 



Let the very natural supposition be maau that tiiu layer 



of liquid contiguous to tlie surfsce of the cilinucr uioves 



with the surface . •Txien the concentric lajers will have c 



motion uur to iriterral friction. There will bo no motiom 



parallel to tlie axis of the; cylinaer, consequently w = C: tne 



only forc^ actin.-i: noi':nil to tiie cylinacir is that cue to li.e 



centrifu^^al forc^j I ut since ti:c velocity of rotation is sriiali 



the prouucts ana squares of x.he velccit\ na:/ be nejlecteu 



without sunsible terror conseqaently u= C also. Thet*efore to 



complettrly aeteraiine the rr.otion of the liquia we need consia- 



cr only the second of the equations for cylinuriccl coorui- 



k- 
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nates. 



Since the initial conaiticn wes one of rest tlie cor.ui- 



tions of equilibriuiii aemand that 

c)v/^t = "^Ir a^/dr^^ 3 dv/ar 4 

r 



2du/lr^d6) di sap J earing since •a= C, 



^ ..2. 



r dv/dz'^aisappearing *ince the motion 



is two uimentionai, 



dv/lr a* ) aisappearing since the motion 



is s:>m'3trical with respect to the axis of tlio cjlinuer, 



Since the an.i^uiar velocity of relation is smell 



2> y/d t bfjc one s a v/ at ; 



ana as soon as tiio iiioticn b-coi.icjs st-^aay 



av/at = C 



Sc that v'ilien the ar.orular velocity is const i*nt ti.c motion of 



the liquia is aeto-rminea from 



r av/ tir^->- b dv/ar = C. 



Solving 



^ 
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rv = Kr-hC/r. 



At an infinite aisience r v 33 finiie consequunllj 



K = r, 



r V = C/r. 



Let 6u bw the an^Ji^SiV velocity of ti;t; cylinuer, ihen 



C = a^lA?/ 



r V = a CO/ r^ 



wrhere a is the raaius of the cylinder. 



But if we consia^r the crlindt^r as filled v^^ith the li- 



quia and consiacr tnu motion of the inclosea liquid 



C = 0, 



auj = K a, 



K = ^, 



r V = r ^^ 



Kence the inclosea liquid moves as a solid boay. 



ISoe Eassei Vol.:.: yV^G^ 23k.) 



j\:ay 2;:,18o'7. 
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